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We review the analysis of the quark and gluon substructures of the scalar mesons from QCD spectral sum rules 
and some low-energy theorems applied to the scalar QCD anomaly current. The present data favour equal 
components of uu + dd and of gg in the wave functions of the low-mass (below 1 GeV) scalar mesons, which make 
the wide a and the narrow /o(980) as 77'dike particles, which can have strong couplings to meson pairs through 
OZI violations. A coherent picture of the other 7 = scalar mesons spectra within this mixing scheme is shortly 
discussed. We also expect the ao(980) to be the lowest isovector ud state, and the Kg (1430) its ds partner. 



1. INTRODUCTION 

Since the discovery of QCD, it has been empha- 
sized JIJ that exotic mesons beyond the stan- 
dard octet, exist as a consequence of the non- 
perturbative aspects of quantum chromodynam- 
ics (QCD). Among these states, the rj meson is a 
peculiar pseudoscalar state, as it should contain 
a large gluon component (the 1st well-established 
gluonium!) in its wave function for explaining 
why its mass is not degenerate with the one of the 
pion but only vanishes in the large N c limit (solu- 
tion of the so-called U(1)a problem) j2), though 
the rj' likes to couple to qq mesons in its decay 
(e.g. rf — > rfKn,...). This peculiar feature has lead 
to certain confusion in the interpretation of its 
nature. Even, at present, many physicists claim 
misleadingly that the rf is a qq state, from the 
conclusion based on the successful quark model 
prediction of its hadronic and 27 couplings. This 
feature is not too surprising due to the manifes- 
tation of the strong OZI-violation for low mass 
state in this channel. However, the rj' mass is 
lower than the direct calculation of the pseu- 
doscalar glueball mass from QCD spectral sum 
rules (QSSR) and lattice QCD (see Table 1). This 
apparent discrepancy can be understood from the 
crucial role played by the rf in the evaluation 
of the U(1)a topological susceptiblity and of its 
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slope p). In the effective Lagrangian approach, 
it has been expected since the pioneer work of 
Nambu-Lonasinio that an analogue of the pion of 
the non-linear a model (NLerM) exists in the lin- 
ear a model (LaM) (J], but this LaM has not 
attracted too much attention in the past because 
it has been observed that the SU(2)l x SU(2)r 
chiral symmetry of QCD is non-linearly realized 
(this has made the success of the chiral perturba- 
tion theory (ChPT) approach Q). It also occurs 
that the Lagrangian of the LaM is not unique 
such that no definite predictions from QCD first 
principles can be done. In addition, the inclusion 
of the resonances into the effective Lagrangian is 
not completly settled ||. 

However, this analogy between the pion and the 
a meson may not be very appropriate as they are 
particles associated to symmetries of a different 
nature (SU(2) L X SU(2) R for the pion and U(l) v 
for the a). Instead, a comparison of the a with 
the 77' meson looks more valuable within a chiral 
U (3) a x U (3)v Lagrangian ^M, as both particles 
are associated respectively to the U(l) axial and 
vector symmetries, where in terms of the quark 
and gluon fields, the corresponding QCD currents 
are respectively, the £^(1)^ anomaly (divergence 
of the singlet axial current): 

M M (s) = tr G afi G ali 

+ "M(475)« > (1) 

u,d,s 
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and the scalar anomaly dilaton current: 

0£ = \/3(a s )G a0 G af3 

+ (l + 7m(a s )) £ m< ^ ' 

u,d,s 



(2) 



which is the trace of the energy- momentum tensor 
6> M „. The sum over colour is understood; q, 
and G^v are respectively the quark, the gluon 
field strength and its dual; m q is the light quark 
running mass; (3 and j m are respectively the (3- 
function and quark mass anomalous dimension. 
In this talk, we shall discuss the a meson from this 
point of view of the analogy with the ?/-meson 
by using QCD spectral sum rules (QSSR) a la 
SVZ H (for a review, see e.g.: and some low- 
energy theorems based on Ward identities. The 
discussions are based on the works in ]Io|-|l2|]. 

2. AN OUTLINE OF QCD SPECTRAL 
SUM RULES (QSSR) 

2.1. Description of the method 

Since its discovery in 79, QSSR has proved to be a 
powerful method for understanding the hadronic 
properties in terms of the fundamental QCD pa- 
rameters such as the QCD coupling a s , the (run- 
ning) quark masses and the quark and/or gluon 
QCD vacuum condensates. The description of the 
method has been often discussed in the literature, 
where a pedagogical introduction can be, for in- 
stance, found in the book In practice (like 
also the lattice), one starts the analysis from the 
two-point correlator: 



Ae" 1 (0\TJ H (x) (J„(0)) f |0) ,(3) 



built from the hadronic local currents Jh{x), 
which select some specific quantum numbers. 
However, unlike the lattice which evaluates the 
correlator in the Minkowski space-time, one ex- 
ploits, in the sum rule approaches, the analyt- 
icity property of the correlator which obeys the 
well-known Kallen-Lehmann dispersion relation: 



dt 



1 



il_ 



t 
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2 Discussions of the qq and/or £ 
states can also be found in e.g. | 
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where ... represent subtraction points, which are 
polynomials in the g 2 -variable. In this way, the 
sum rule expresses in a clear way the duality 
between the integral involving the spectral func- 
tion IvmpH (t) (which can be measured experimen- 
tally), and the full correlator ipH(q 2 ). The latter 
can be calculated directly in the QCD Euclidean 
space-time using perturbation theory (provided 



that 



-<r 



(m being the quark mass) is 



much greater than A 2 ), and the Wilson expan- 
sion in terms of the increasing dimensions of the 
quark and/or gluon condensates which simulate 
the non-perturbative effects of QCD. 

2.2. Beyond the usual SVZ expansion 

Using the Operator Product Expansion (OPE) 
||, the two-point correlator reads: 

Mq 2 ) - £ ' 



3=0,2,4,... I y I 

C{q\v){0{v)) , 



(5) 
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where v is an arbitrary scale that separates 
the long- and short-distance dynamics; C are 
the Wilson coefficients calculable in perturbative 
QCD by means of Feynman diagrams techniques; 
{0{v)) are the quark and/or gluon condensates 
of dimension D. In this paper, we work in the 
masslcss quark limit. Then, one may expect the 
absence of the terms of dimension 2 due to gauge 
invariance. However, it has been emphasized re- 
cently |m| that the resummation of the large or- 
der terms of the perturbative series, and the ef- 
fects of the higher dimension condensates due e.g. 
to instantons, can be mimiced by the effect of 
a tachyonic gluon mass, which might be under- 
stood from the short distance linear part of the 
QCD potential. The strength of this short dis- 
tance mass has been estimated from the e + e~ 
data to be ||,|§: ^A 2 ~ -(0.06 - 0.07) GeV 2 , 
which leads to the value of the square of the 
(short distance) string tension: a ~ — |a s A 2 ~ 
[(400±20) MeV] 2 in an (unexpected) good agree- 
ment with the lattice result |^7j of about [(440 ± 
38) MeV] 2 . The strengths of the vacuum con- 
densates having dimensions D < 6 are also un- 
der good control: namely 2m(qq) = —m%f% from 
pion PCAC, (a s G 2 ) = (0.07 ± 0.01) GeV 2 from 
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e + e — > / = 1 data |25| and from the heavy quark 
mass-splittings @, a s (g<?) 2 ~ 5.8 x 1CT 4 GeV 6 
||, and 5 3 (G 3 )w 1.2 GeV 2 (a s G 2 ) from dilute 
gaz instantons |13J . 

2.3. Spectral function 

In the absence of complete data, the spectral 
function is often parametrized using the "naive" 
duality ansatz: 

~Tmi/> H (t) ~ 2M]Pffj5(t - Mfj) + 

7T 

"QCD continuum" x 6(t — t c ), 

(6) 

which has been tested |J using e + e~, r-decay 
data, to give a good description of the spectral 
integral in the sum rule analysis; fjj (analogue to 
f„) is the the hadron's coupling to the current ; 
2n is the dimension of the correlator, while 
is the QCD continuum's threshold. 

2.4. Form of the sum rules and optimiza- 
tion procedure 

Among the different sum rules discussed in the lit- 
erature within QCD || (Finite Energy Sum rule 
(FESR) §§, r-like sum rules @,...), we shall 
mainly be concerned here with: 

• The exponential Laplace unsubtracted sum rule 
(USR) and its ratio: 

f°° 1 
CJt) = / dt t n exp(-tr) - lmtp H (t) , 
Jo t 

n„ = -— log£ n , (»>0); (7) 
dr 

• The subtracted sum rule (SSR): 

[°° dt 1 
£_i(r) = — exp(-ir) - Imip H (t) 

Jo t 7T 

+^ff(0) ■ (8) 

The advantage of the Laplace sum rules with re- 
spect to the previous dispersion relation is the 
presence of the exponential weight factor, which 
enhances the contribution of the lowest resonance 
and low-energy region accessible experimentally. 
For the QCD side, this procedure has eliminated 
the ambiguity carried by subtraction constants, 
arbitrary polynomial in q 2 , and has improved the 
convergence of the OPE by the presence of the 



factorial dumping factor for each condensates of 
given dimensions. The ratio of the sum rules is 
a useful quantity to work with, in the determina- 
tion of the resonance mass, as it is equal to the 
meson mass squared, in the usual duality ansatz 
parametrization. As one can notice, there are "a 
priori" two free external parameters (r, t c ) in the 
analysis. The optimized result will be (in princi- 
ple) insensitive to their variations. In some cases, 
the ^-stability is not reached due to the too naive 
parametrization of the spectral function. One can 
cither fixed the t c - values by the help of FESR (lo- 
cal duality) or improve the parametrization of the 
spectral function by introducing threshold effects 
fixed by chiral perturbation theory, in order to 
restore the instability of the results. The results 
discussed below satisfy these stability criteria. 

3. UNMIXED GLUONIA MASSES AND 
DECAY CONSTANTS 

Before discussing the specific scalar channel, let's 
present the situation of gluonia/glueball mass cal- 
culations as a guide for the forthcoming discus- 
sions. 

3.1. The currents 

In addition to the pseudoscalar and scalar cur- 
rents introduced previously, we shall deal with 
the tensor and 3-gluon currents (standard nota- 
tions): 

J3 = dfabcGapGp^G^a > (^) 

3.2. Masses and decay constants 

The unmixed gluonia masses from the unsub- 
tracted QCD Spectral Sum Rules (USR) @ are 
compared in Table 1 with the ones from the 
lattice J27],|3ll in the quenched approximation, 
where we use the conservative guessed estimate 
of about 15% for the different lattice system- 
atic errors (separation of the lowest ground states 
from the radial excitations, which are expected 
to be nearby as indicated by the sum rule anal- 
ysis; discretisation; quenched approximation,...). 
One can notice an excellent agreement between 
the USR and the lattice results, with the mass hi- 
erarchy: M ++ < Mq-+ M 2 ++, expected from 
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Table 1 



Unmixed gluonia masses and decay constants from QSSR |12| compared with the lattice. 





Name 


Mass [GcV] 






Vt~c [GeV] 


Sa [MeV] 






Estimate 


Upper Bound 


Lattice [27,31 






0++ 


G 


1.5 ±0.2 


2.16 ±0.22 


1.60 ± 0.16 


2.1 


390 ± 145 




3G 


3.1 


3.7 




3.4 


62 


2++ 


T 


2.0 ±0.1 


2.7 ±0.4 


2.26 ±0.22 


2.2 


80 ± 14 


o-+ 


P 


2.05 ±0.19 


2.34 ±0.42 


2.19 ±0.32 


2.2 


8-17 



some QCD inequalities However, this is not 
the whole story ! Indeed, one can notice that in 
the pseudoscalar channel, the predicted value of 
the mass of the v is too high compared with 
the mass of the rj' , which is not surprising for the 
reasons explained in the introduction. We shall 
see in the next section that the same phenomena 
occur for the scalar channel. 

4. UNMIXED SCALAR GLUONIA 

4.1. The need for a low mass <jb from the 
sum rules 

Using the mass and decay constant of the scalar 
gluonium G in Table 1 from the USR_[Eq||), into 
the SSR (Eq.|) sum rules, wher e fl3| ^ s (0) ~ 



T6(/3i/7r)(a s G 2 ), one can notice |K,|l2j that one 
needs a second resonance f\ with a lower mass 
the as) for a consistency of the two sum rules 
' Using, e.g., M„ B ~ 1 GeV |, one gets @,||: 
fa B as 1 GeV, which is larger than fa — -4 GeV. 

4.2. Low-energy theorems (LET) for the 
couplings to meson pairs 

In order to estimate the couplings of the gluonium 
to meson pairs, we use some sets of low-energy 
theorems (LET) based on Ward identities for the 



3 All previous old sum ru e analysis based on a one- 
resonance parametrization 3^ l£ad_to inconsistencies, as 
emphasized many years ago in |l 0^1 2[ ] . 

4 The introduction of the 1/q 2 in the OPE has been also 
shown to restore the mass hierarchy between the gluonium 
and meson scales p6) , and then improve further the con- 
clusion reached here. A recent analysis using gaussian sum 
rules and instanton approach also favours ^two-resonance 
parametrization of the spectral function fejfl . 

5 We cannot fix both the mass and decay constant of the a 
from the two sum rules. We give in the original papers a 
more complete analysis for different values of the a mass 
ranging from 500 to 1 GeV. 



vertex: 

V(q 2 ^(p-p') 2 =0) 
V'(0) 



ee (H(p)mH(p')) 



2m 
1 , 



H > 



(10) 



and write the vertex in a dispersive form. H can 
be a Goldstone boson (ir, K,r)%), a r\\- U(1)a- 
singlet , or a as- Then, one obtains the sum 
rules for the hadronic couplings: 



1 E 

4 ^ 



9SHH 



V2fs 



2M 2 H , 



X gsH H V2f s /M 2 s 



4 

<r B ,<r' ,G 



(11) 



• Decays into 7T7t: Neglecting, to a first approxi- 
mation the G-contribution, the <jb and a' B widths 
to 7T7T, KK,... (we take M a > as 1.37 GeV as an 
illustration) are Q fj 



T(a B 



7T7r) 
7T7r) 



0.8 GeV 
2 GeV , 



(12) 



which suggests a huge OZI violation and seri- 
ously questions the validity of the lattice results in 
the quenched approximation. Similar conclusions 
have been reached in [p~5| |l8[| . For testing the 
above result, one should evaluate on the lattice, 
the decay mixing 3-point function U(0) responsi- 
ble for such decays using dynamical fermions. 
• Decays into rj'rj and rjm Using rj' as cosBpr\\ 
and r\ as sinOprji, where 9p is the pseudoscalar 

"However, one can notice that the a coupling to pion pairs 
decreases like i.e., a too low value of the mass say 

less than 500 MeV, leads to a width less than 100 MeV. 
This feature does not favour a too low value of the a mass 
and questions the interpretation of some data. 
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mixing angle, the previous LET implies the char- 
acteristic gluonium decay (we use Mq w 1.5 
GeV and assume a G-dominance in the sum rule) 

CoE3i: 



r(G 
r(G 



7777 ) 



r(G 



(5 - 10) MeV , 
0.22 : 



7777') 

g Gm ~ sin 6» P gGr/ri' ■ (13) 

• Decay into 47r°: Assuming that the G de- 
cay into 47T° occurs through <7_b0\b, and using the 
data for /o(1.37) — > (47r°)g, one obtains from the 
previous sum rule ]To|Jl^| : 

r(G -> ciijcri} -> 4tt) w (60 - 140) MeV . (14) 

4.3. 77 widths and J/-0 — > 7^ radiative de- 
cays 

These widths can be estimated from the quark 
box or anomaly diagrams 0,^2| ■ The 77 widths 
of the a, a' and G are much smaller (factor 2 to 5) 
than T(r/ — > 77) ~ 4 keV, while B(J/tfj — > 7 <r, cr' 
and G) is about 10 times smaller than B(J/ip — > 
777') w 4 10 -3 . These are typical values of gluonia 
widths and production rates juj . The absence of 
the a in 77 scattering and its presence in J/\& 
radiative decays Q are a strong indication of its 
large gluon component. 

5. UNMIXED SCALAR QUARKONIA 
5.1. The a (980) 

The ao(980) is the most natural and economical 
meson candidate associated to the divergence of 
the vector current: d ll V^{x) = {rn u — rnd)u(i^^)d. 
Previous different sum rule analysis of the as- 
sociated two-point correlator gives 0]: M aa ~ 
1 GeV and the conservative range f ao ~ (0.5 — 
1.6) MeV (f n = 93 MeV), in agreement with 
the value 1.8 MeV from a hadronic kaon tadpole 
mass difference approach plus a ao dominance of 
the KK form factor. A 3-point function sum rule 
analysis gives the widths fllT 21 12]]: 

r(a -> r/ir) ~ 37 MeV , 

r(a ^77) ~ (0.3- 1.5) keV , (15) 
while from SU (3) symmetry, we expect to have: 



9a K+K 



The value of this coupling to KK reproduces 
present data jl7],[l8|]. Analogous sum rule anal- 
ysis in the four-quark scheme ||l],|9| gives simi- 
lar values of the masses (see also the lattice re- 
sults in p2j| ) and hadronic couplings but implies 
a too small value (compared with the data) of 
the 77 width due to the standard QCD n 2 
loop-diagram factor suppressions. Such a factor 
have not been taken properly in the literature. 
The {uu — dd) quark assignement for the do (980) 
is supported by present data and alternative ap- 
proaches p7l-pl. 



5.2. The isoscalar partner S2 = uu + dd of 

the d (980) 

Analogous analysis of the corresponding 2-point 
correlator gives Mg 2 « M ao as expected from a 
good SU (2) symmetry, while using 3-point func- 
tion and SU (3) relation, its widths are estimated 

to be lilJril: 



T{S 2 ^n+TT-) 
T(5 2 



120 MeV , 



77) - y r ( fl o ' 
~ 0.7 keV 



77) 



(17) 



5.3. The ^(1430) = ds and S 3 = ss states 

An analysis of the K^-ao mass shift due to SU(3) 
breakings (strange quark mass and condensate) 
131 fixes the mass of the KX to be around 1430 



MeV. If a candidate around 900 MeV is con- 
firmed, it will then be hard to reconcile with a 
qq structure. An analysis of the S3 over the Kq 
2-point functions gives 112]: 



M S JM K5 
M S3 
fs 3 



1.03 ±0.02 = 
1474 MeV , 
(43 ± 19) MeV 



in agreement with the lattice result 
3-point function leads to |L2j: 



T(S 3 -> K+K-) 

r(s 3 -> 77) 



(73 ± 27) MeV 
0.4 keV . 



(18) 
while the 



(19) 



la.o'tf'R • 



(16) 



In the usual sum rule approach (absence of large 
violations of the OPE at the sum rule stability 
points), one expects a small mixing between the 
S*2 and 5*3 mesons before the mixing with the glu- 
onium ob- 



G 



5.4. Radial excitations 

The propreties of the radial excitations cannot be 
obtained accurately from the sum rule approach, 
as they are part of the QCD continuum which 
effects are minimized in the analysis. However, 
as a crude approximation and using the sum rule 
results from the well-known channels (p,...), one 
may expect that the value of \ft~ c can localize ap- 
proximately the position of the first radial excita- 
tions. Using this result and some standard phe- 
nomcnological arguments on the estimate of the 
couplings, one may expect fll2||: 



M S' 2 

T(S' 2 - TT+O 

nS' 2 - 77) 
nS' 3 - 77) 



1.3 GeV, 

(300 ± 150) MeV, 
(4 ± 2) keV , 
1.7 GeV, 
(112 ±50) MeV, 
(1 ± .5) keV. 



(20) 



5.5. We conclude that: 

• Unmixed scalar quarkonia ground states are not 
wide, which excludes the interpretation of the low 
mass broad a for being an ordinary qq state. 

• There can be many states in the region around 
1.3 GeV (a 1 , S3, S' 2 ), which should mix non- 
trivially in order to give the observed /o(1.37) 
and /o(1.5) states (see next sections). 

• The /j(1.7) seen to decay mainly into KK |l8| , 
if it is confirmed to be a ++ state, can be the 
first radial excitation of the S3 = ss state, but 
definitely not the pure gluonium advocated in |23| . 

6. SCALAR MIXING-OLOGY 

Many scenarios have been proposed in the liter- 
ature for trying to interpret this region P|,|l7|]~ 
p3| . However, one first needs to clarify and to 
confirm the data (H| for a much better selection 
of these different interpretations. 

6.1. Mixing below 1 GeV and the nature 

of the a and Jo (980) 
In so doing, we consider the two-component mix- 
ing scheme of the bare states (ct_b, S'2): 



l/o >= 

1(7 > = 



- smt 
cos( 



Wb > 



cost 
sint 



\S 2 > 

\s 2 > 



A sum rule analysis of the off-diagonal 2-point 
correlator 1 35 |Jl^] : 



i J d 4 x e lqx (0\T/3(a s )G 2 (x) 
x m M (0)|0), 



(22) 



responsible for the mass-shift of the mixed states 
gives a small mass mixing angle of about 15°, 
which has been confirmed by lattice calculations 
using different input for the masses |^3|] and from 
the low-energy theorems based on Ward identi- 
ties of broken scale invariance Jl6|, if one uses 
there the new input values ]9|,^5| of the quark and 
gluon condensates. In order to have more com- 
plete discussions on the gluon content of the dif- 
ferent states, one should also determine the decay 
mixing angle. In so doing, we use the predictions 
for gb, Si — * 77 obtained in the previous sec- 
tions and the data r(/ — > 77) ~ 0.3 keV. Then, 
we deduce a maximal decay mixing angle and the 
widths HQ: 

~ (40-45)°, 
< 134 MeV , 
« 2, 

w (300 - 700) MeV , 
w (0.2 - 0.5) keV . (23) 



r(/ ^+o 

9f K+K- /ff/oTT + TT- 

r((T — > 7T + 7T~) 

r(cr -» 77) 



The huge coupling of the /o to KK or ss comes 
from the gluon component through the large mix- 
ing with the a. For this reason, the fo can have a 
large singlet component, which is also suggested 
by independent analysis |T^,|l5|. Extending the 
previous J /if) — > 7+ X analysis into the case of the 
6, one obtains the new result within this scheme 



Br[cj) 7 + / (980)] w 1.3 x 10" 



(24) 



(21) 



in good agreement with the Novosibirsk data of 
(1.93 ± 0.46 ± 0.5) x 10~ 4 . In the same way, the 
large coupling of the /o(980) to KK could also 
explain its production from D s — » 37r decay [ ^7| . 

6.2. Mixing above 1 GeV and nature of the 

/ (1.37), / (1.5), / (1.6) and / (1.7) 
As already mentioned previously, this region is 
quite complicated due to the proliferation of 
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states. We shall give below some selection rules 
which can already eliminate some of the different 
schemes proposed in the literature: 

• The / (1.37): If it decays into aa -> (47r°) s , it 
signals mixings with the a, a' and G. 

• The / (1.5): 

— If it decays into aa and 7/77, this signals a gluon 
component. 

— If it also couples to tttt and KK, this signals a 
qq component which may come from the S' 2 , S3. 
Then, it can result from the qq mixings with the 
a, a' and G, like the / (1.37). 

— If it couples weakly to tttt and KK, while 
the ratio of its r\r\ and 7/77 is proportionnal to 
l/sin 2 #p, then it can be an almost pure gluo- 
nium state, which can be identified with the G 
in Table 1 obtained in the quenched approxima- 
tion (this approximation is expected to be much 
better at higher energies using 1/N C arguments 

EH)- „ 

• The Jo (1.7): If it decays mainly into KK p[, 
it is likely the radial excitation S' 3 of the Ss(ss) 
state rather than a gluonium. 

7. CONCLUSIONS 

We have seen that there are increasing experimen- 
tal evidences for the existence of the a and other 
scalar mesons. QCD spectral sum rule (QSSR) 
and low-energy theorem (LET) provide a first 
analysis of their gluon and quark substructure 
beyond the usual effective Lagrangian approach. 
Present data favour a maximal quarkonium- 
gluonium (qq-gg) mixing scheme for the a and 
/o(980) mesons. Analogous scenarios, though 
more complicated, are encountered for higher 
mass scalar states. Within the scheme, the wide 
a and the narrow /o(980) appear to be ?/-like par- 
ticles, which can have strong couplings to meson 
pairs through OZI violations. The a (980) is a 
ud state, while its strange partner ds is expected 
to be the ATq(1430) of PDG. More experimen- 
tal tests are needed for selecting different phe- 
nomenological schemes. We wish that in the near 
future our understanding of the scalar mesons will 
be improved further. 
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